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Abstract 

We consider the effect of spatial correlations on sources of polarized electromagnetic radiation. The 
sources, assumed to be monochromatic, are constructed out of dipoles aligned along a line such that their 
orientation is correlated with their position. In one representative example, the dipole orientations are 
prescribed by a generalized form of the standard von Mises distribution for angular variables such that 
the azimuthal angle of dipoles is correlated with their position. In another example the tip of the dipole 
vector traces a helix around the symmetry axis of the source, thereby modelling the DNA molecule. We 
study the polarization properties of the radiation emitted from such sources in the radiation zone. For 
certain ranges of the parameters we find a rather striking angular dependence of polarization. This may 
find useful applications in certain biological systems as well as in astrophysical sources. 

1 Motivation 

In a series of interesting papers Wolf 0, ||, || studied the spectrum of light from spatially correlated sources 
and found, remarkably, that in general the spectrum does not remain invariant under propagation even 
through vacuum. The phenomenon was later confirmed experimentally [§, [5| ^| and has been a subject of 
considerable interest fil]. Further investigations of the source correlation effects have been done in the time 
domain theoretically |8| and experimentally [^| . Several applications of the effect have also been proposed 
]l0| , 11, |l2|, 13, [I!]]. In a related development it has been pointed out that spectral changes also arise due 
to static scattering [jl5], [l6|, [l7], 18, [l9| and dynamic scattering [^0[ ^l], |2^, |2^, 24 . In the present paper we 
investigate polarization properties of spatially correlated sources. Just as we expect spectral shifts for such 
sources, we expect nontrivial polarization effects if the correlated source emits polarized light. 

The basic idea can be illustrated by considering two polarized point sources Pi and P 2 which are located 
along the z-axis at a distance 2z apart, in analogy to a similar situation considered by Wolf Q for the 
unpolarized case. The electric field at the point Q located at large distances Ri and R2 from these point 
sources can be written as, 

B = e(fi,w)— — + e(r 2 ,u)— — (1) 

Here r% and r*2 are the position vectors of the points Pi and P2 . We calculate the coherency matrix for this 
electric field at the point Q, 



< ETE* > = 



< e*(r 1 ,uj)e j (f 1 , ui) > < e|(f 2 , w)ej(f 2 , w) > 



R\ " R\ 



R1R2 



h.c. (2) 



We are interested in investigating the effect of the third term which arises due to cross correlation between 
the sources P\ and P2 ■ It is clear that this term will give nontrivial contribution to the polarization in the far 
zone. In the present paper we investigate the contribution of this term by considering only monochromatic 
waves. It turns out to be useful to understand this simple idealization before treating the realistic case of 
partially polarized radiation. 

In order to illustrate the contribution of the cross correlation term we consider the situation where the 
two sources are simple dipoles, pi and p 2 located at z and — z respectively and are oriented such that their 

1 current address: Physics Department, University of Rochester, Rochester, NY 
2 e-mail: pkjain@iitk.ac.in 



1 



polar angles 9x = 62 = 9 P and the azimuthal angles 0i — ~02 = tt/2. We assume that 9 P lies between 
and 7r/2. The strength of the dipolcs is po and they radiate at frequency u>. We compute the electric field 
at point Q located at coordinates (i?, #,</>), as shown in Fig. 1, such that R » z. The electric field in the 
far zone is obtained by the addition of two vectors pi ■ RR — p\ and pi ■ RR — pi with phase difference of 
2uiz cos9/c. The vector p ■ RR — p at any point (R, 9, 4>) is ofcourse simply the projection of the polarization 
vector p on the plane perpendicular to R at that point. We keep only the leading order terms in z/R in 
computing the total electric field. 

The observed polarization is obtained by calculating the coherency matrix, given by 

j-csfi sf*sn (3) 



v (e;e,) {E t E tl , 

The state of polarization can be uniquely specified by the Stokes's parameters or equivalently the Poincare 
sphere variables p8| . The Stoke's parameters are obtained in terms of the elements of the coherency matrix 
as: So = J x\ + ^22) Si = Jn — J22, S2 = J12 + J21, S3 = i(Jzx — Jxi)- The parameter So is proportional 
to the intensity of the beam. The Poincare sphere is charted by the angular variables 2%, and 2ip, which can 
be expressed as: 

Si = So cos 2\ cos 2-0, S2 = So cos 2\ sin 2-0, S3 — So sin 2% (4) 

The angle \ (~ 7r /4 < X < 7r /4) measures of the ellipticity of the state of polarization and tp (0 < ip < tt) 
measures alignment of the linear polarization. For example, \ = represents pure linear polarization and 
X — 7r/4 pure right circular polarization. 

The Stokes parameters at the observation point Q for the two dipole source are given by 

S a = [4 cos 2 (wz cos 6»/c) cos 2 6» p sin 2 6* + 4 sin 2 (wz cos 6»/c) sin 2 6» p (cos 2 6» sin 2 + cos 2 0)] 

Sx = ( [4 cos 2 {loz cos 9/c) cos 2 9 p sin 2 9 + 4 sin 2 (luz cos 9/c) sin 2 9 P (cos 2 9 sin 2 - cos 2 0)] 



.2, 
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' ' ' 8 sin 2 (u>z cos 9/c) sin 2 9 P cos 9 sin <fi cos < 



c 2 R 

S3 = ( . ° ] 8 cos(ujz cos 9/c) sin([jz cos 9/c) cos 9 P sin 9 P sin 6* cos 
V c R J 

The resulting polarization in the far zone is quite interesting. At sin0 = 0, 9 = n/2 the wave is linearly 
polarized (x = 0) with ijj — 0. As 9 decreases from ir/2 to 0, x > and the wave has general elliptical 
polarization with = 0. At a certain value of the polar angle 9 = 9t the wave is purely right circularly 
polarized. As 9 crosses 9t, the linearly polarized component jumps from to tt/2, i.e. 2ip changes from to 
7r. The value of the polar angle 9t at which the transition occurs is determined by 



tan(wz cos 9 t /c) = ± sin 9 t / tan 9 p (5) 



From this equation we see that as z — > 0, the transition angle 9 t is close to zero for a wide range of values of 
9 p . Only when 9 p — > 7r/2, a solution with 9 t significantly different from can be found. In general, however, 
we can find a solution with any value of 9t by appropriately adjusting z and 9 P . For sin0 > 0, we find 
qualitatively similar results, except that now the linear polarization angle 2ip increases smoothly from as 
9 goes from tt/2 to 0. Furthermore the state of polarization never becomes purely circular, i.e. the angle 2x 
is never equal to tt/2 although it does achieve a maxima at 9 close to the transition angle 9t given by Eq. 5. 

We therefore find that for wide range of parameters the polarization in the far field region shows a 
nontrivial dependence on angular position which is governed by the cross correlation term. In general this 
term will give nontrivial contribution unless the two sources have identical polarizations. 
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2 A One Dimensional Spatially Correlated Source 

A simple continuous model of a monochromatic spatially correlated source can be constructed by arranging 
a series of dipoles along a line with their orientations correlated with the position of the source. The dipoles 
will be taken to be aligned along the z axis and distributed as a gaussian exp[— z 2 /2a 2 ]. The orientation of 
the dipole is characterized by the polar coordinates 9 p ,<f> p , which are also assumed to be correlated with the 
position z. A simple correlated ansatz is given by 

exp [a cos(6>p) + (3z sin(0 p )] 

Ni{a)N 2 {pz) 1 ' 

where a and (3 are parameters, N\{a) — t:Iq{q) and N2((3z) — 2ttIq((3z) are normalization factors and 
Iq is the Bessel function. The basic distribution function exp(acos(# — 9q)) used in the above ansatz is 
the well known von Mises distribution which for circular data is in many ways the analoque of Gaussian 



distribution for linear data |25L E6L 27 . For a > this function peaks at 9 = 8q- Making a Taylor 
expansion close to its peak we find a gaussian distribution to leading power in 9 — 9 . The maximum 
likelihood estimators for the mean angle 9o and the width parameter a are given by < sin(6> — 9$) >= and 
< cos(# — #o) >— dlog(Io(a)) / da respectively. In prescribing the ansatz given in Eq. 6 we have assumed 
that the polar angle 9 p of the dipole orientation is uncorrelated with z and the distribution is peaked either 
at 9 P = (ir) for a > (< 0). The azimuthal angle 4> p is correlated with z such that for (3 > and z > 0(< 0) 
the distribution peaks at cf> p = ir/2(3ir/2). 

We next calculate the electric field at very large distance from such a correlated source. The observation 
point Q is located at the position (R, 9, </>) (Fig 1), measured in terms of the spherical polar coordinates, and 
we assume that the spatial extent of the source a « R. The electric field from such a correlated source at 
large distances is given by, 



/OO pTT pI'K 

dzexp(~z 2 /2a 2 ) / d9 p / d<t> p 
oo JO JO 
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cxp (a cos 9 P + (3z sin <j> p ) 
2^I {a)I {(3z) 



cxp (ibjzil ■ z/cj x (p ■ RR — p) (7) 



where p is a unit vector parallel to the dipole axis, po is the strength of the dipole, lu is the frequency of 
light and Iq denotes the Bessel function. Since we are interested in the radiation zone we have dropped all 
terms higher order in z/R. The resulting field is ofcourse transverse i.e. E ■ R = 0. We have also assumed 
that all the dipoles radiate at same frequency and are in phase. The correlation of the source with position 
is measured by the parameter (3. 

It is convenient to define scaled variable z = z/a, X = X/a where A = 2ituj/c is the wavelength, and 
(3 = (3a. The integrations over 9 P and (f) p can be performed analytically. We numerically integrate over z 
for various values of position of the observation point, the parameter a which determines the width of the 
distribution of 6 p and for different value of the correlation parameter f3. 

We first study the situation where f3 > and a > 0. The result for several values of (9, <j>) are given in 
figures 2,3 which show plots of the Poincare sphere variables 2% and 2ij). The scaled wavelength X ~ X/a 
of the emitted radiation is taken to be equal to ir, i.e. the effective size of the source a is of the order of 
the wavelength A. The results show several interesting aspects. The ellipticity of the state of polarization 
shows significant dependence on the position of the observer. The angle \ = 0, i.e. the beam has pure linear 
polarization, for the polar angle cos(0) = 0, 1 for all values of azimuthal angle (f>. It deviates significantly 
from as cos(0) varies from to 1. For sin(^) = 0, 2% = ir/2 at some critical value 9t as cos(#) varies 
between and ir/2, i.e. the state of polarization is purely right circular at 6 = 9 t . For sin(</>) > 0, 2\ also 
deviates significantly from and displays a peak at some value of 9. The precise position of the peak is 
determined by the values of the correlation parameters a and (3. 

The alignment of linear polarization also shows some very interesting aspects. For sin(</>) = 0, we find 
that V is either or 7r depending on the value of 9. The transition occurs at the same critical value of 9 
where the angle \ shows a peak. The state of polarization is purely linear with the electric field along the 9 
for cos(#) = and then acquires a circular component for increasing values of cos(0). At the transition point 
9 = 9 t , the polarization is purely circular. With further increase in value of 9 the state of polarization is 
elliptical with the linearly polarized component aligned along <j>. The transition point is clearly determined 
by the condition Si = Jn — J22 = 0. 

For other values of sin((/>) we find ip — for cos 6* = and then deviates significantly from as 9 
approached 9 t , finally levelling off as cos 9 approches 1. The final value of -0 at cos 9 = 1 depends on the 
correlation parameters and sin^ but for a wide range of parameters 2ip > ir/2. We see from Fig. 3 that the 
linear polarizations from sources of this type shows striking characteristic, i.e. that the polarization angle ip 
is either close to or tt/2 depending on the angle at which it is viewed. 

For s'mcf) < the Poincare sphere polar angle 2\ is same as for sin</> > 0, however the orientation of the 
linear polarization 2ip lies between w and 27r, i.e. in the third and fourth quadrants of the equatorial plane 
on the Poincare sphere. For a particular value of <f> the azimuthal angle ip{<j)) = —ip(—(j)). 

If we change the sign of a we do not find any change in linear polarization angle ip however the value of 
X changes sign, i.e. the state of polarization changes from right elliptical to left elliptical. Change in sign of 
(3 also leaves tp unchanged while changing the sign of x- Changing the signs of both a and (3 produces no 
change at all. 

In the case of the limiting situation where (3 = we find, as expected, linear polarization is independent 
of the angular position, i.e. x = and ip = 0. This is true for any value of the parameter a, which determines 
the polar distribution of the dipole orientations. Hence we see that the effect disappears if either the effective 
size of the source a = or the correlation parameter (3 = 0. The effect also dissappears in the limit a — > 00. 
In this limit the distribution of 9 p is simply a delta function peaked at and hence our model reduces to a 
series of dipoles aligned along the z-axis, which cannot give rise to any nontrivial structure. In the numerical 
calculations above we have taken the effective size of the source a of the order of the wavelength A. If the 
size a « A, the effect is again negligible since the phase factor ivzR ■ z/c in Eq. 7 is much smaller than one 
in this case. 

Hence we find that in order to obtain a nontrivial angular dependence of the state of polarization the 
size of the source, assumed to be coherent, has to be of the order of or larger than the wavelength as well as 
the correlation length 1/(3. 
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2.1 Transition angle 



From our results we see that there exists a critical value of the polar angle 9 at which the state of linear 
polarization changes very rapidly. This is particularly true if we set sin0 = where we find that that the 
orientation of linear polarization tp suddenly jumps from (or tt) to n/2 at some critical value of the polar 
angle 9 = 9 t . We study this case in a little more detail. The 9 and <f) components of the total electric field 
is given by, 

c z R Io{oe) 

£ = j^PO c -i"(t-R/c) 2sinha A 

c 2 R airlo(a) 

A = ^ dze- z2 / 2 °\m(ouzcos9/c)t^- 
J-oo h{pz) 

In this case the Stokes parameter S2 = 0. For (3 > (<)0, S3 > (<)0 and hence \ > (<)0- The point where 
the polarization angle 2ip jumps from to 7r is determined by the condition S\ = 0. This is clearly also the 
point where 2\ — ±n/2. Explicitly the condition to determine the critical value 9 t is, 

A 2 = <j\ 3 a 2 e-° 2 " 2c ° s2e </ c2 sin 2 9 t h{a) * . 

2sinh a 

This can be used to determine 9 t as a function of a, (3. The result for cos6> t as a function of a is plotted 
in figure 4 for several different values of (3. For any fixed value of the parameter (3, the transition angle 9 t 
decreases from tt/2 to as a goes from zero to infinity. This is expected since as a becomes large the polar 
angle distribution of the dipole orientations, peaked along the z axis, becomes very narrow and hence the 
resultant electric field is aligned along the z axis for a large range of polar angle 9. Furthermore we find, as 
expected, that as (3 goes to zero the transition angle also tends towards 0. 



3 Helical Model 

We next study an interesting generalization of the model discussed above. Instead of the having the peak 
of the <j)p distribution fixed to — n/2 for z < and tt/2 for z > we allow it to rotate in a helix circling 
around the z-axis. In this case we replace the <j) p dependence by exp[/3(</f) p — £z)]. As z goes from negative 
to positive values, the peak of the distribution rotates clockwise around the z-axis forming an helix. This is 
a reasonable model of the structure of DNA molecule and hence has direct physical applications. We study 
this in detail by fixing the azimuthal angle of the dipole orientation <p p = £z and the polar angle 9 p to some 
constant value, i.e. the <p p and 9 P distributions are both assumed to be delta functions. This allows us to 
perform the z integration in Eq. 7 analytically. The resulting state of polarization, described by Poincare 
sphere angles 2\ and 2-0 are shown in Figs. 5-8 . In this model we can extract a simple rule to determine 
the transition angle for the special case 9 P = tt/2 and £ = nir where n is an integer. We set sin0 = for 
this calculation since it is only for this value that the polarization becomes purely circular for some value of 
9 = 9 t and the linearly polarized component flips by tt/2 at this point. A straightforward calculation shows 
that this transition angle 9 t is given by: 

cos 2 9 t = n\/2 

Here n represents the number of ir radians that are traversed by the tip of electric field vector along the 
helical path and A is the wavelength. In order to get at least one transition A < 2/n. In the special case 
under consideration there is atmost one transition. However in general the situation is more complicated 
and for certain values of 9 P and £, more than one transitions are possible. Some representative examples are 
shown in Figs. 5-8. 
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4 Conclusions 



In this paper we have considered spatially correlated monochromatic sources. We find that at large distance 
the polarization of the wave shows dramatic dependence on the angular position of the observer. For certain 
set of parameters the linearly polarized component shows a sudden jump by ir/2. If the symmetry axis 
of the source is taken to be the z-axis, the polarization shows a sudden transition from being parallel 
to perpendicular to the symmetry axis of the source, as the polar angle is changed from ir/2 to 0. The 
sources considered in this paper are idealized since we have assumed coherence over the entire source. 
For small enough sources, such as the DNA molecule, this may a reasonable approximation. In the case of 
macroscopic sources, this assumption is in general not applicable. However in certain situations some aspects 
of the behavior described in this paper may survive even for these cases. For example, we may consider a 
macroscopic source consisting of large number of structures of the type considered in this paper. As long as 
there is some correlation between the orientation of these structures over large distances we expect that some 
aspects of the angular dependence of the polarization of the small structures will survive, even if there does 
not exist any coherent phase relationship over large distances. Hence the ideas discussed in this paper may 
also find interesting applications to macroscopic and astrophysical sources. As an interesting example we 
consider astrophysical sources of radio waves. It is well known that the polarization angle of these sources is 
predominantly observed to be aligned either parallel or perpendicular to the source orientation axis |p9|| . This 
difference has generally been attributed to the existence of different physical mechanism for the generation of 
radio waves in these sources. Our study, however, indicates that this difference in observed polarization angle 
could also arise simply due to different angles of observation. Hence orientation effects must be considered 
before attributing different physical mechanisms for differences in observed polarizations of these sources. 
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sin((|>) = 
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sin((|)) = .75 

^ 2 I 1 1 1 




Figure 2: The polar angle on the Poincare sphere 2%, which is a measure of the eccentricity of the ellipse 
traced by the electric field vector. For pure linear polarization 2\ — and for pure right circular polarization 
2% = 7r/2. The 3-D plot shows 2\ as a function of cos# and sin-0 where 8 and <j> are the polar and azimuthal 
angles of the point of observation. The 2-D plots on the right show the corresponding slices of the 3-D plots 
for different values of sin0. The upper and lower plots correspond to /3 = 1, a = 0.25 and (3 = a = 1 
respectively. 
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Figure 3: The azimuthal angle on the Poincarc sphere 2ip. This measures the orientation of the linearly 
polarized component of the wave. The 3-D plot shows 2ip as a function of cos 6* and sin^ where 9 and 
(j) are the polar and azimuthal angles of the point of observation. The 2-D plots on the right show the 
corresponding slices of the 3-D plots for different values of sin <j>. The upper and lower plots correspond to 
P = 1, a = 0.25 and (3 = a = 1 respectively. 
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Figure 4: The critical value of the polar angle at which the state of linear polarization shows a sudden 
transition for sin <j> = as a function of the parameters a and (3 which specify the distribution of the dipolc 
orientations. For any given value of the parameters a and (3, electric field is parallel (ip = 0) to z axis if the 
cosine of the observation polar angle cos 9 is less than cos 9 t . On the other hand electric field is perpendicular 
to the z axis if cos 9 is greater than cos t . 
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Figure 5: The polar angle on the Poincare sphere 2% (radians) for the helical model as a function of cos(#) 
(A = 0.2n,9 p = tt/2,£ = tt). 
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Figure 6: The azimuthal angle on the Poincare sphere 2ip (radians) for the helical model as a function of 
cos(0) (A = 0.2tt, 9 P = it/2, £ = tt). 
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Figure 7: The polar angle on the Poincare sphere 2% (radians) for the helical model as a function of cos(0) 
(A = 0An,e p = tt/4, £ = tt). 
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Figure 8: The azimuthal angle on the Poincare sphere 2ip (radians) for the helical model as a function of 
cos(0) (A = 0.4tt, p = tt/4, £ = tt). 
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